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ABSTRACT 

In this report the main physics processes in LH2 tank during prepress and rocket flight are 
studied. The goal of this investigation is to analyze possible hazards and to make risk 
assessment in proposed LH2 tank designs for SLS with 5 engines (the situation with 4 engines 
is less critical). For analysis we use the multinode model (MNM) developed by us and presented 
in a separate report and also 3D ANSYS simulations. We carry out simulation and theoretical 
analysis the physics processes such as (i) accumulation of bubbles in LH2 during replenish 
stage and their collapsing in the liquid during the prepress; (ii) condensation-evaporation at the 
liquid-vapor interface and tank wall, (iv) heating the liquid near the interface and wall due to 
condensation and environment heat, (v) injection of hot He during prepress and of hot GH2 
during flight, (vi) mixing and cooling of the injected gases due to heat transfer between the 
gases, liquid and the tank wall. We analyze the effects of these physical processes on the 
thermo- and fluid gas dynamics in the ullage and on the stratification of temperature in the liquid 
and assess the associated hazards. A special emphasize is put on the scaling predictions for 
the larger SLS LH2 tank. 
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Nomenclature 


p 

= density 

P 

= dynamical viscosity 

p 

= pressure 

Pr 

= Prandtl number, Pr=pCp/k 

Pc 

= critical pressure for gas 

Ra 

= Raleigh number 

T 

= temperature 

Re x 

= Reynolds number 

T c 

= critical temperature for gas 

Gr 

= Grashof number 

T s 

= surface or interface temperature 

Nu 

= Nusselt number 

u 

= velocity 



c 

= sound velocity 



M 

= Mach number, M = u/c, M 0 = u/c 0 



C v 

= specific heat for the constant volume 



C P 

= specific heat for the constant 

Subscripts: 

y 

pressure 

9 

= gas 

9 

= ration of specific heats; y= Cp/c v 

L 

= liquid 

R 

= gravity 

env 

= environment 

R, 

= radius 

t 

= tank 

s, 

= radius of tank 

w 

= tank wall 

Ssu 

= surface of tank cross-section 

ull 

= ullage 

r 

= total surface of tank 

eva 

= evaporation 

Ld 

= thermodiffusion length 

P 

=condensation 

9l 

= heat of evaporation 

con 


h T 

= heat transfer coefficient 

d 


j 

= mass flow density [j = up] 



J 

= mass flow [ J = upS] 



q 

= heat flow density 

Abbreviations: 

Q 

= heat flow 



Qc 

= convection heat flow 

LH2 

= Liquid Hydrogen 

Qr 

= radiation heat flow 

GH2 

= Gaseous Hydrogen 

Pt 

= thermal expansion coefficient 

MNM 

= Multinode Model 

5 

= boundary layer thickness 

SLS 

= Space Launch System 

k 

= thermal conductivity 

3D 

= Three dimensional 

V 

= kinematic viscosity 




Introduction 


2 


We study main physical processes in the LH2 tank during the prepress and rocket flight. The 
goal of this research is to analyze the hazards associated with the functioning of the external 
cryogenic LH2 tank of the Space Shuttle (Fig.1) during the pre-pressurization process and in the 
course of the flight and their scalability to the proposed SLS design with 5 engines (the situation 
with 4 engines is less critical). We use both the protocol, material parameters, and the telemetric 
data for Shuttle flights [1, 2] to fit an uncertain parameter of the multinode model (MNM), 
developed in our group with use results of the works [3-5]. After an appropriate rescaling MNM 
is further applied to study the temperature stratification, the ullage pressure dynamics and the 
associated risks in the LH2 tank in the framework of the proposed SLS designs. 



H,»29 4m 


Fig.l. Design of the Shuttle tanks: the shape and parameters of the LH2 tank 


1. Boiling near the tank wall and bubbles’ collapsing: scaling of the 

bubbles’ volume 

LH2 is characterized by a low temperature of boiling [6-8]. Therefore, a very weak heating of 
the liquid boundary layer near the tank wall by the environment heat flow Q gnv =q e nvSsur leads to 

boiling LH2 (Fig. 1-1) by creating bubbles. We show below that the total volume of the bubbles is 

2 

proportional to the tank height squared H ( due to increase both of the power Q env and the rising 
time of the bubbles. 
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Parameters: 

*Ht=29.4mand 
I S=774m 2 ; 

J R,= 4 m; 

1 P t =2nR f : 

I p L =70.37kg/m i : 
q L =4.5x10-J/kg; 

1 G=9660J/kg/K; 

1 KpO.I W/K/m; 

I p,=1.1x10' 5 Pa.s; 
I T,=20.4K; 
f P r = 0.023K' 1 ; 

I R H )=4124J/kg/K; 
l Pv=p/(Rh2Tl); 

^ p=1atm. 
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Fig.l-l. Boiling near the LH2 tank wall due to the environmental heat flux. 


1 .1 . Nucleate boiling of liquid in the boundary layer near the wall 


During the replenish stage bubbles are generated near the tank wall due to the heat flux q 
from the hot environment. The bubbles mass generation rate can be assessed as 

_ 


m 


bubble 


<//. 


, ™ bubble ~ 0.3 5frg / sec at q env = 200w / m 


(D 


The velocity of bubbles rising under the buoyancy force is given by 

= eV. Jo, - o )- ^ 

m 


Pvapoi^bub ^ S^hubi^PL P vapor) ^ P hub ’ 




( 2 ) 


bub 


v huh 


3 C, 


' 0(0.4 + 1.8)-, (0.04 <C,< 0.4) 


Here we estimated the typical size of the bubble is by the LH2 capillary 
radius r capi „ = j2a LH1 ! gp L = 2 mm (cr IJi2 = \.l •XQr’N I m) . The bubbles number density along 
the tank wall can be approximated using 


dP bubble _ dp, 


bubble / bubble ,, staf 

V bub 


_ qJsui 


' Phuhble ( x ) “ 


Qw^xur 

slat 

9 1 y bubble 


-x 


dt dx q, 

Eq. (3) implies 

y _ f PbuhhJJ 1 — yO 

V bubble ~ | r bubble 

0 Phubble . 0 

It follows that the bubble volume is a quadratic function of the tank height length H t . 


'JL* 


(3) 


(4) 


1 . 2 . 


Collapsing of the bubbles during the prepress 


During the initial stage of the prepress bubbles will collapse due to increase of the pressure 
from 1 amt to 3.2atm over ~16 sec. The characteristic time-scale of the collapsing can be 
estimate from the energy balance 


(T - T ) 

^ r L K L~ - Nu ,r b ■' conaps = KubPs^L’ 


'bub 


where p s vapor is density of saturated H2 vapor at p=3atm and 


(5) 


Nusselt number Nu /r = 0.332 

• r hn 


.tfa/ 

Pi. ^ bubble ' bubble, 0 




J 


C_lPl 


vl/3 


K 


50, 


L J 


The bubble surface temperature is determined by the ullage pressure [9] 


We find from Eq. (5) 


t. 


collapse 


T>{p) = T c 


y /53 

p_ 

\Pcj 


, T s = 25 K at p- 3.2 atm 


P vapotP bubble, 0 . 

3 K l (T a -T l )Nu’ 


<lsec at p>\.5atm 


(6) 


(7) 


Thus, the bubbles will collapse completely during of the initial stage of the He injection, 
leading to increase of the ullage volume. The ullage and bubble volumes in Space Shuttle are 
V bU bbie=58.25ft 3 and Vuii a ge = 372 ft 3 (see Tablel), i.e. V bub bie / V u n ag e =16%. While, in SLS LH2 tank 
with 5 engines V bubb | e = 162.25ft 3 and V u n ag e=620 ft 3 , i.e. V bubb | 6 / V u n age =26%. This effect has to be 
taken into account in the simulation of LH2 tank and in particular of the prepress regime. 


2. Main processes of heat transfer in the LH2 tank 

The general heat balance in the ullage LH2 tank (Fig. 2-1) can be written as 

Qiryect Qen\’ Qullage + QwaU + Qi ■ ( 8 ) 

Here Q mject +Q em is the total heat added to the ullage tank as a result of hot gas injection 

from the diffuser and the interaction with environment. This heat is spent on heating the gas 
inthe ullage, 
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V„ ullage volume 
p ullage pressure 

Diffusor T ullage temperature 



Fig.2-1. LH2 tank design. Main physical processes and variables 


Q u „ , heating the tank wall, Q wa „, and on condensation-evaporation processes including 

heating of the surface layer of the liquid , Q, . The heat flows associated with these processes 

are sketched in Fig. 2-1, which also defines the variables of the physical model. 

We now proceed to analyzing these pathways and estimating their relative importance for 
the problem of temperature stratification and pressure dynamics in the ullage. 


2.1. Condensation-evaporation at the liquid-vapor interface 


A possible risk associated with the condensation-evaporation processes at the liquid-gas 
interface is a deviation of the ullage pressure from its permissible bounds (Fig.1). Here we 
analyze the evaporation of LH2 due to the hot He or GH2 injection from the diffuser. We notice 
that the change of the ullage pressure is slow (~sec) compared to the timescale of 
thermalization of the liquid-vapor interface (< 0.1msec). As a consequence, the quasi- 
equilibrium near the interface is established for all times with the surface temperature T s 


determined by the vapor pressure p vapour = p GH2 [9] : 


f 


> 1/A 


T = T. 


Poo , 

Pc 


Ps(T) = p c 


j.'N 


(9) 


Here T c =33. 2K and p c =13.2atm are the critical temperature and pressure of hydrogen, p s is 
the pressure of saturated vapour near the interface and X=5.3 for H2 (see Fig. 2-2) [6-8], 
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Fig.2-2 Dependence of saturated vapor pressure on the surface temperature of LH2 


Note, that p s =1atm at T S =T L =20.4K. The partial H2 vapor pressure is p GH2 a 3.2atm=47psi and 
p G H2 ls 2.25atm=32psi during the prepress and flight, respectively. For such pressures the liquid 
surface temperature T s > 71. Since T >T, and T>T the convection of both liquid and vapor 

is absent. (Moreover, our 3D simulations have also shown that the convection induced by the 
hot gas injection from the diffuser and natural convection near the wall is relatively small near 
the gas-liquid interface.) Thus, the heat balance in condensation-evaporation processes is 
established mainly by means of thermo-diffusion at the given surface temperature T s (p GH 2 ) . The 
thermodiffusion is driven by the equation 

dT _ d 2 T 
L - gP '- 8 dt ~ Kl * dt 2 ’ 

The solution of is 

£r 

L * pL * dt ~ Kig dt 2 ’ 


T(0,t) = T s = const, T(r,co) = T l0 (10) 

7X0,0 = T s = const, T(r, oo) = T L0 


T(r,t) = T L0+ {T 1 -T L0 )erf 


4k L t 


■\ 


J 


Here p, g ,C, i ,,rc l K are the density, heat capacity, and thermo-conductivity of liquid (gas- 

vapor); x=0 is the liquid-vapor interface, erf is the error function. Using Eq. 
Error! Reference source not found, we find that the heat flux from the heated surface to the 
interface 


il L (x,t) = -K L 


0T(x,t ) 


dx 


= I CrpMr, -r,.) ! 

V 7Ct 


T 


T-T, 


‘ L L 


'DL 


( 11 ) 


where L DI = 


- HElL; 


I c ,.Pl 


is thermodiffusion lenght and Q, = S, cj, . 


This makes it possible to estimate the condensation mass flow j cond from the following 
balance equation: 


a,,u = (q L + <\ jnz (r s -r)) = Ql - Qs = S L " t - — 

V \ l ) 

,fK^c7o, ( 

jconj(‘) = 


-T 


yjK L C,p L (T s — 7l)— \JK G h iPgh 2 {jg ~ T s ) 


(0 


Dg\ ) y 


(12) 


Ql + C G hi {j g T) 


1 



Here the first term in the left part of Eq. (12) describes the condensation heat flow rate and 
the second describes the evaporation flow rate. The mass of condensed vapor is given by (1 1) 

/ 

M cond {t)=\S L {t)j cond {t)dt (13) 

o 


According to Eq. (12) condensation will be dominant whenever 

(T,-UPoo,))<(T,(Po a .)-T,), 


K I ^ I Pi 


^GOX C G()x P G Q } 


’> TSPgo,) ~ T c 


Y /5 - 3 

Pc, Ox 

Pc ) 


(14) 


According to the condition (14), the condensation is dominant when the gas temperature 
near the liquid-vapor interface T g < 70K during the prepress when p==3.2 atm and T g < 65K 
during the flight when p=2.25atm. For the estimation of the thermo-diffusion length in the gas 
we used characteristics of pure H2 prevalent in the vicinity of the interface. This result agrees 
with the experimental data [10] presented in Fig. 2-4. The point in Fig. 2-4, right correspond of 
our estimation for p=32psi. 

MNM simulation showed that the effective temperature of the ullage gas near the interface 
is less than 70K. Therefore, the condensation is dominant in the LH2 tank during the prepress 
(Fig. 2-5) when the partial GH2 pressure is higher than 2atm (Fig. 2-6). However, it is expected 
(see estimation in Fig. 2-4, left) that the condensation mass flux is relatively small in comparison 
with rates of injection of both hot He (j He =0.6kg/sec) and H2 (j H2 =0.47kg/sec) from the tank 
diffuser. We note that the condensation rate found in MNM simulation is an overestimate since 
a higher averaged temperature in a relatively wide layer is used to calculate the rate. 



Fig.2-3: Temperature profile near the liquid-gas interface 



Fig. 2-4 Evaporated mass (left) and critical temperature (right) as function of vapor He pressure: 
theoretical estimations (left at p G n2 = 2.25atm) and experimental data [9] 
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Fig. 2-5. Time-traces of the condensation flow rate (left) and the condensed mass (right). 

Results of the MNM simulation. 



Fig. 2-6. Partial pressures of H2 vapor and He near the liquid-vapor interface. Results 

of the MNM simulation. 


2.2. Cooling of the ullage gas by condensation at the interface 

The heat flux into the liquid associated with the condensation is given by 

(T-T s ) 


d(J. 


dt 


V') 


■S L (t) = (T-T s )S L (t) 


k £iRa 


nt 


u Dg 




is thermodiffusion lenght in the gas 


The total heat due to the condensation of GH2 can be estimated as 


QA 0 * (T g - T s )S L Jrc L C L p L (4t/x), 


(15) 


(16) 


QgL ( ^prepress ) * 4.6 X 10 V and QJt^) *10 V 

The heat of injected He during of the prepress (t=100sec) and the total heat of injected hot 
GH2 during the flight (t fl ig h t =500sec) are estimated as follows 


Qinjecl ~ C p.He^He^ He 1 prepres ~ 10 J, Qi„ iec , ~ C p GH 2 Tq H 2 J CH it flight 9 . 1-10 J . 


inject p.GH 2 1 GH V* GH 2 * fl ighl 


(17) 
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Thus, the condensation weakly affects the general ullage heat 

balance: QgjM flight ) ^ Qinject — 1% The effect of condensation strengthens considerably due 
to the sloshing (see Section 5). 


2.3. Evaporation near the wall-liquid-gas interface 


The hot injected GHe or GH2 heats the tank wall causing heat flow along the wall. This flow 
will induce evaporation of liquid near the wall-liquid-gas interface (Fig 2-7). The heat flow is 
Given by the following expression: 


Qwl = 2xR,d u 


dT 
— — 

ax 


\ 


J /interface 


*2 vR,d w KAT s ) 


(T ws -T s ) 


Ad.. 


(18) 


The evaporation flow J and the evaporated mass M near the wall-liquid-gas 

evap.w evap.w 

interface can be estimated from the heat balance equation 


, -2*Rd K I . (^-^) , f _r J 

dgvap.w^L evap.w ~~ . , > ^evap.w ~ J 

^*ws Ql ^*ws 


dt 




Interface 

•80 -60 40 -20 0 20 4G N 90 

Tmtts 


Heat flow 
along wall 

100 


MNM simulation 



T, +3K 


after prepress start 


^Interface 


Y inch 


(19) 


Fig. 2-7 Heat and mass evaporation flows near the wall-liquid-gas interface (left); Wall temperature 
stratification (results of 3D, in center, and of MNM, right, simulation). 


It follows from the MNM and 3 D simulation 


dT\ W 

K w — <4000— [(T -T )< 30 K, 

^ ^/interface ™ ™ 


Kws^lTW/m/K and Ad > 0.15m], Therefore, according to Eq. (19) 


•4«p.h < 0. 8g / sec, w < 0.04kg during the prepress (20) 

Note that the total mass of gas in the ullage is greater than 16kg. The condensed vapor mass 
accumulated during the prepress is ~ 0.4kg (see Sec. 2.1). Thus the effect under consideration 
is negligibly small during the flight due to >T (Fig. 2-7, right). 


2.4. Heating of the tank wall by the hot gas: the main mechanism of the 

gas cooling 


The heat flow from the hot gas in the ullage at T g to the tank wall at T w can be written as 


10 



( 21 ) 


* *=- = M r « - r»), V = ^ )' 


: 25 




1+ 


f 0.492 V" 6 


Pr 


EF 

/w 2 /:’ 

\ -J6/9 


Here we use the correlation formulas for the turbulent regime (see Section 4.2.1) as the 
Rayleigh number Ra>3x 10 10 and Re x >10 5 for x>1m. The heat transfer coefficient „ and the 

temperature length l =K,/h do not depend 


'gw 1 

on 


( k * 4.5 mm for p = 2Aatm, T g = 300AT, T w = 120 K). 
The total heat flux is 


- T ^ SS u,(0 = <C w (T w )p w SS u f 

df <? k -T w ) C w (T w )p w d w 

dt r^TJ ’ ^ w) ^ 


(22) 


The temperatures of the ullage and the tank wall depend on the coordinate x along the wall. 
For the numerical simulations the temperature dependences of the specific heat C w and thermal 
conductivity for Al-Li 2195 alloys have been used. The changes in the temperature of ullage gas 
and in the temperature of the tank wall during the prepress are shown in Fig. 2-8. The solid lines 
are plots of MNM-simulated time-traces of temperatures near the top (upper red curve) and 
bottom of the tank (blue). The wall temperature TJt) increases quickly during the initial stage 

due to the small value of C w (T L )=18J/kgK at T=23K that increases strongly to 840J/kgK at 
T=255K. 




Fig. 2-8 Time traces of the ullage and wall temperature during the prepress. 

Heating of the tank wall is the main mechanism of the hot gas cooling in the ullage. 


3. Ullage pressure and temperature stratification 


3.1 . Fitting of the MNM and He pulses parameters to the Space Shuttle 

data 


1 1 


The values of main parameters at the end of replenish (EOR) stage and at the engine start 
command (EST) are shown in Table 1 (from [1]). The variation of ullage pressure during the 
period, t = -103sec to t = 100, is shown in Fig. 3-1 [1]. The pre-pressurization starts at 

^plpress ~ -102.8sec with injection of the hot He and finishes at t=-86.5sec. The nominal He flow 

rate is J He =1-3lb/sec =0.6kg/sec and He temperature is T= 572.4R=318K. The nominal 
parameters of supply of gaseous He or H2 for STS -115 are presented in Table 1. 



Table 1 : Nominal LH2 tank (STS 115) 


EOR 


Ambient Pressure: 

1.024atm 

14.65 psia t 

Ullage Pressure: 


0.20 psig 

Ullage Pressure: 


14.85 psia 

Liquid Level: 


1044.6X1 

Ullage Volume: 

10.5m 3 

372 ft 3 

Ullage Mass: 


27 lb m 

Total Tank Volume: 

1498m 3 

52.915 ft 3 

Bubble Volume: 

1.65m 3 

58.25 ft® 

Liquid Volume: 

1486m 3 

52.484 ft 3 

Vented LH2 Density: 


4.41253 lb m 

LH2 Propellant Load: 


231.588 lb m 


Maximum He flow J He =1 .3lb/sec=0.6kg/sec 


Maximum H2 flow J H2 =3x0.35lb/sec=0.476kg/sec 


ESC Through T-0 

LH2 Losses Prior to ESC: 

1 13 lb m 

Ullage Pressure: 

32.35 psig 

Ullage Pressure: 3.24atm 

47 psia 

GHe (Anti-Ice) Mass Added: 

1 lb m 

GHe (Prepress) Mass Added: 

38 lb m 

Ullage Volume @ ESC 27.1Sm 3 959 ft 3 

Ullage Mass (a ESC 

68 lb m 

ET Liquid Volume @ ESC 

52.320 ft 3 

ET Liquid Mass @ ESC 

231.475 \b m 

Press’d LH2 Density 

4.42427 lb m ft 

Predicted LH2 Usage to T-0: 

1.792 lb m 

LH2 Transferred to SSMEs: 

>8 lb* 

Loaded LH2 at T-0: 

229.625 lb n , 

Predicted Ullage Volume - °o of 

1.83 

Liquid Volume 

Required Ullage Volume - °o 

1.54 

of Liquid Volume 


The ullage pressure drops at t>t ( ^ press due to cooling of the hot injected He (see Sec. 2.3) 

and expansion of the ullage volume under the pressure increase. Pulses of He injection are 
applied when the ullage pressure drops below the threshold value of p th =20.54 psig=35.2psia. 
The form of the He mass flow pulses is shown in Fig. 3-2. In more detail, the pulses have 
amplitude j He =0.3lb/sec=0. 136kg/sec and their form is given by: 

[ 1 -exp(-(r -/,.)/ T fmn ), t<t, + r 

jV)=j„e\ , ' (23) 

[ exp(-(/ - f, - T p ) / T lail ) , t>t'+T p 

where at f=f, the ullage pressure drops to the threshold value of p, h and i-pulse is generated. 
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Fig. 3-2: The form of the He mass flow pulses. 


Ullage temperature stratification and gas velocity distribution during the first stage of the 
prepress found in 3D simulation are shown in Figs. 3-3 and 3-4. Temperature distributions 
during the prepress and first stage of the flight are shown in Fig. 3-5. One can see that 
temperature smoothly decreases from the high temperature near the tank top (T=300K) to the 
low temperature at the interface (T=23K) during the prepress, and the temperature distributions 
becomes more pronounced during the flight (Fig. 3-5). 

We chose a single fitting parameter of MNM, numerical coefficients (3 U in the velocity 
correlation for the turbulent natural convection (see Eg. (25) ), in such a way as to fit the 
temperature distribution found in the 3D simulation (Figs. 3-3 - 3-5). The sawtooth modulations 
of the ullage pressure induced by the He pulses were fitted to the telemetric data. For the fitting 
we chose parameters of He impulses in Eg. (23) (The pulse duration r and its front and tail 

pulse 

durations t and t ). 

front tail 


h? = O.l50„K L 
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8, = 0.565x[(l + 0.494 Pr 2/3 ) / Ra x Pr 23/I5 ] oc jc 7/ '° 

Amplitude of He mass flow was 1 .3lb/sec =0.6kg/sec. The fitting gave the following values 
for the parameters: 

• p u = 0.2 

• t = 0.12sec, t = 0.4sec, and t = 0.4sec 

front tail pulse 
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Boundary Layers: 60 layers inflation , the first layer 0.1mm , rate 1.05 


Fig. 3-3: Distribution of the ullage temperature (upper) and velocity (below) during the 
first stage of the prepress: 3D simulation at t=10sec (left)-and 16sec (right) after the 

prepress start. 




Fig. 3-4: Temperature stratification during the prepress: Results of the 3D simulation for 
different times after the beginning of prepress. 
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Fig. 3-5: Temperature stratification during the prepress and the first stage of the 
Shuttle flight: Results of 3D simulation for different times after prepress start. 

Stratification of the ullage and wall temperature found in a MNM simulation is shown in Fig. 
3-6. Results of fitting time-traces of the ullage pressure and temperature near the tank top, 
where the thermo-sensor is located, are shown in Fig. 3-7. We see that the results of the MNM 
simulations fit well the STS-1 15 data. We emphasize that the only fitting parameter of the MNM 
was the numerical coefficient |3 U in the correlation relation for the natural convection velocity (Eq. 
(24)). This coefficient depends on the fluid’s properties and is essentially independent on the 
tank design. Therefore we used MNM with this fitting parameter for analysis of the rescaled LH2 
SLS tank. 



Fig. 3-6. Stratification of the ullage (left) and wall temperature (right) during the pre- 
press and the initial stage of the flight. Results of the MNM simulation. 
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Fig. 3-7: Time-traces of ullage pressure (left) and ullage temperature (right) near the 
tank top during the pre-press and the initial stage of the flight. Results of MNM 
simulation. 


3.2. Scaling considerations for the LH2 tank for SLS with 5 engines 


The MNM with fitting parameter (B u = 0.2 are further used in modeling the prepress and fight 
regimes for the proposed design of the SLS LH2 tanks to assess the associated risks. Results 
of these simulations for 5/3 scaling SLS tank are shown in Figs. 3-8 and 3-9. Fig. 10 illustrates 
the ullage temperature distribution in 5/3 scalable SLS Lh2 tank found from 3D simulation for 
prepress stage. We increased the mass flows of injected He and GH2, and also ullage volume 
and tank volume, due to its extension, in 5/3 times. 

It can see from comparison of Figs. 3-4 and 3-6 with Figs. 3-8 and 3-10 that the temperature 
stratification is more pronounced in the 5/3-rescaled SLS LH2 tank due to a lesser mixing of the 
initial cold and injected hot gases. Results presented in Figs. 3-7 and 3-9 indicate that the time 
traces of the ullage pressure in Space Shuttle tank and in 5/3-rescaled LH2 tank are close. 
However the numbers of the sawtooth variations of the ullage pressure initiated by injection of 
He impulses during the prepress are different because of difference in the GHe masses and 
surfaces of the tank wall corresponding to the ullage volumes in Space Shuttle and SLS tanks. 




Time (s) 


Fig. 3-8. Stratification of the ullage (left) and wall temperature (right) during the pre-press 
and the initial stage of the flight for 5/3 scalable SLS Lh2 tank. Results of MNM 

simulation. 
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Fig. 3-9: Time-traces of ullage pressure (left) and ullage temperature (right) near the 
tank top during the pre-press and the initial stage of the flight for 5/3 scalable SLS LH2 

tank. Results of the MNM simulation. 

We see that the observed time trace of the ullage pressure for SLS 5/3-rescaled LH2 tank falls 
within the permitted bounds (see Fig. 3-11) 
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Fig. 3-10: Temperature stratification in 5/3 scalable SLS Lh2 tank: Results of 3D 
simulation for different times after the beginning of the prepress. 
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Fig. 3-11. Requirements to the ullage pressure [1] 
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4. Temperature stratification in the liquid 

The temperature stratification patters in the tank is shaped by two main physical processes: (i) 
condensation-evaporation on the gas-liquid interface and (ii) natural convection along the tank 
wall (Figure 4-1). 



1. Condensation-evaporation - 


T 5 =23.8K 


GH2 

p GH2 =2-25atm 

T e =50-300K 
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2. Turbulent natural convection: 
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Fig.4-1. Physical processes contributing to the temperature stratification in cryogenic 

tank 


4.1 Heating of liquid in the boundary layer near the gas-liquid interface 

The liquid layer near the liquid-gas interface is heated due to condensation. The surface 
liquid temperature is determined by the partial pressure of H2 vapor p GH2 near the interface (Eq. 
(9)). It can be seen from Fig. 2-6 that the partial GHe pressure near the interface 
p G H2 < 32psi=2.25atm. Therefore, the surface liquid temperature, according to Eq. (9), is T s < 
23. 8K (T s -T l < 3.4K). The width of this overheated liquid layer is determined by the 
thermodiffusion length (Fig. 4-1) of the liquid (Eq. (11)) which increases with time and 
approaches a value of L dl = 1.5 cm toward the end of the flight. This conclusion explains results 

of 3D simulation presented in Fig. 4-2 that shows temperature stratification in the liquid near the 
interface in the end of the flight (tfiig ht =500sec). 

The upper bound for the mass of the overheated liquid in this layer can be estimated as 

while the total LH2 mass in thefilled tank is 105,000kg. 


The heat flow into the liquid associated with the condensation is given by 


Q l - dQl =K, 

dt L L DL (t ) 


^ r ^S L (t) = (T s -T L )S L (t\ 


Kt 


The total heat associated with the condensation of GH2 is 


&(') * (T, - ' x). &('«.) “ 5 x 10 5 J 

while the total heat of GH2 injected during the flight (t fl jg ht =500sec) is 


(25) 


(26) 
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(27) 


Q ullage ~ ^^Ox^GOX^ ullage ~ ^ xOx 
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tank 
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= 2.5-10 V 


Q m]ea = = 3.4 • 10 V, J„„ = 1 3kg / sec, 7\ = 800/? = 444/: 


" /jCbr* inj 17 inf flight 


mj 


Thus, the additional heat of this overheated layer has a relative small value: 

a(w) / a f »= ai% - 


4.2 Heating of liquid in the boundary layer near the wall 


4.2.1. Natural convection correlations 

There is a constant heat flow across the tank wall q w ~200W/m 2 . The natural convection flow 
generated by this heat is highly turbulent in a large tank such as LH2 tank. With characteristic 
Reynold number Re~10 6 . As a consequence a theory dealing with the natural convection in the 
tank is a theory of turbulent flows. Due to the notorious closure problem (Figure 4-2) little 
progress can be done by a rigorous mathematical analysis of the equations of motion and 
additional physical assumptions have to be made which cannot by derived from the equations of 
motion and the boundary conditions. These assumptions have to be validated by comparison of 
the theoretical predictions with the experiment. We note that the results of this subsection are 
applicable to natural convection in both the liquid and the ullage space, provided the 
Boussinesq approximation (Figure 4-2) hold. 
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The problem of turbulent natural convection 


Boussinesq approximation: 

p 0 (l + fiT')«. p 0 T=T-T m 

d x u 2 4 - dyTU — vdpJ 4 - ilgT' = 0 . 



P = 

Averaged momentum 
tra nsfer eq-n 
( Navi er-S takes): 

Averaged energy 
Transfer eq-n: 

Continuity eq-n: 


0 x uT 4 - d y vT = xtfTT 

Oj-U + dyV = 0. 


Boundary conditions: v w = 0. = 0. Toe, and or T" ir 


- n 

Closure problem: u- ^ « ,uv ^Tiv, etc. 

Fig. 4-2. The closure problem for the turbulent flow 

In contrast to the forced convection there seems to be no consensus in the literature about 
what are the correct assumptions for the natural convection theory [11-16], In what follows we 
outline the physical picture which is correct in our opinion, derive certain correlation relations 
based on this picture and validate the predictions by comparison to experimental evidence and 
high fidelity simulations which we believe to be correct. 

To set up the stage for the discussion, let’s consider Figure 4-3, which depicts a turbulent 
flow near the wall due to the natural convection. One should distinguish two characteristic 
length-scales of this flow: the viscous or conduction length-scale I, and the turbulent length 
scale 5. The former lengh-scale corresponds to the inner sublayer of the flow, where the heat 
and momentum are transferred by thermodiffusion and viscous fiction, respectively (for LH2 the 
Prandtle number is ~1, therefore both the thermal diffusion and the viscous friction are 
characterized by the same length-scale). 
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Turbulent natural convection. Length-scale separation 


Conduction (viscous) length scale: 



/ <£ s 


/ \ 

"'0.5 mm "20cm 

(Shuttle, LH2 tank) 


Fig.4-3. Length-scales separation in the turbulent natural convection 

The latter length-scale, 5, corresponds to the outer region of the turbulent boundary layer, 
where the transport of both heat and momentum is carried out by the turbulent pulsations. For 
well-developed turbulent flow as in the LH2 tank, there is a huge separation between the two 
length-scales with ratio l/5~10' 3 . 

The turbulent heat and momentum transport is highly efficient. Therefore, we expect that the 
all the heat absorbed by the viscous sublayer from the wall is transferred into the outer turbulent 
sublayer (Figure 4-4). This is a very powerful assumption and it leads to very different results 
from those obtained by assuming that rather the total momentum flux from the wall is 
transferred into the outer sublayer, as happens in forced convection and as was assumed in an 
influential paper by Eckert and Jackson [12]. 

Assuming the heat flux q w into the outer region, where the microscopic (as distinguished 
from the turbulent) viscosity and thermo conductivity are negligible, we are left with the following 
physical variables which can determine the characteristic velocity and temperature scales in the 
outer region: heat flux, q w , gravity acceleration g, volumetric expansion coefficient, specific heat 
C L , density p and streamwise coordinate along the wall x. Out of these variables one can build 
just a single velocity scale u and a single temperature scale T' = T — T x (T~ is the temperature 

in the bulk of the liquid far away from the wall) [13]: 


u-A 


n V /3 
C lP , 


T'-B 


f \ in 


j 


{gpx) 


- 1/3 


(28) 


where A and B are empirical constants, which have to be determined from an experiment. 
Remarkably, the temperature is seen to decrease streamwise! This counterintuitive prediction 
can be understood if one realizes that the heat absorbed from the wall is diffused by the 
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turbulent pulsations in ever increasing width of the turbulent outer layer 5 down the stream, 
leading to a smaller temperature increase per same heat absorbed. 



Constant heat flux assumption 





Assumption / <C 6 


On the dimensional grounds [*], the turbulent 
characteristic (“average) velocity Is given by 


u =£ a ( gag* 


Empirical constant 
— (functions of Pr) 


And the temperature (deficit) T*=T-T m 


{gfa) 


-1/3 

Empirical constant 
(functions of Pr) 


The average temperature in the turbulent outer layer decreases streamwlse I 


* W. K. Georgs and S. P. Capp, A theory for natural convection turbulent boundary layers next to heated 
vertical surfaces, Int. J. Heat Mass Transfer 22, p.813 (1979) 


Fig.4-4. Assumption of constant heat flux into the outer turbulent region. 

To calculate the actual x-dependence of 5 we consider the heat balance for a certain control 
volume in the boundary layer (Figure 4-5): 


q w x&C , puT'S 


(29) 


Eqs. (29-30) imply that the width of the outer region grows linearly with the streamwise 
coordinate x, 5~x. This fast growth explains the decrease of the characteristic temperature in 
the turbulent region with the height of the tank x. 
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Estimation of the turbulent boundary layer width 


The heat balance equation: 



qx ~ Cl puT'S 

t \ 

Heat flux into the Heat flux out of 
control volume the control volume 



S = x(AB )- 1 

Tbe turbulent characteristic (^averags) temperature decreases streamwise: the turbulence 
stirs the absorbed heat in the ever wider boundary layer streamwise: ~ x 


Fig. 4-5. Estimation of the turbulent boundary layer width. 

A very important consequence of the predicted streamwise temperature decrease is that the 
main temperature change occurs in the inner conductance region of width l w . In fact, near the 
wall the temperature changes linearly with the distance from the wall y: 


Cj> ' dy 


(30) 


So the change of the temperature over the width l w of the conductance region 

A Th-2^- (31) 

C,pK 

is independent of x, i.e., for sufficiently large x it becomes much larger than the 
characteristic turbulent outer region temperature scale T'~x~ 1/3 . It follows that for large x the 

temperature of the wall is T w « AT + T X , , which sets the characteristic scale of the temperature in 

the conductance region. This is a very important conclusion which means essentially the 
temperature increase of the boundary layer due to the heat flux is independent of the height of 
the tank for sufficiently large tanks (Figure 4-6). 

The foregoing argument can be considered as a theoretical justification of the constant heat 
transfer coefficient for the natural convection h, which is a very well established empirical 
correlation relation, known in the literature [17]: 
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(32) 
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+ (0.492/Pr) 97 * 6 ] 169 


An important contribution of our theoretical analysis is assertion that not only is the wall 
temperature constant for the constant heat flux turbulent natural convection, but that the heating 
essentially occurs in the conductance sublayer only, i.e., that AT = T W -T X . This latter property 

allows one to calculate not only the characteristic temperature in the inner sublayer but its width, 
too, which is important for the stratification analysis. 


.onstant heat transfer coefficient 


For large x the temperature gradient is substantial only in the viscous sublayer of width I. 
This can explain the constant heat transfer coefficient for turbulent boundary layer. 



Empirical constants 


in the inner (viscous/conduction) sublayer : 

q or T w - r x 

% K ^ 

C L p Oy ! 

q = h(T ir — TJx;) 


i.e., the heat coeff icient h is independent of x, 

which explains a well-established empirical heat-transfer 

correlation: 
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M'(Pr) = [l t (0.492 / Prrj 1 " 5 


Fig. 4-6. Constant heat transfer coefficient. 

The precise value of the temperature drop near the wall AT =T W - T„, including the 

numerical factor, can be calculated using the empirical expression for the heat transfer 
coefficient (33). 

To summarize the main predictions of the foregoing theoretical arguments: 

♦ The absorbed heat is accumulated in the turbulent outer sublayer 

♦ The width of the turbulent sublayer increases streamwise 

♦ As a consequence the average temperature in the turbulent sublayer 
decreases streamwise 
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• The temperature gradient is substantial in the conductance (viscous, 
inner) boundary sublayer 

* The temperature is streamwise constant in the conductance boundary 
sublayer 

The important implication of these predictions, if proved (experimentally/numerically) correct, 
for the temperature stratification in the tank is the following: The temperature of the warm layer 
accumulated beneath the gas/liquid interface is independent of the tank height. In other words, 
the temperature stratification pattern is essentially invariant with respect to the tank rescaling. 
By “essential” we mean heating which is likely to be of any consequence from the Net Positive 
Suction Pressure Condition (NPSP) perspective. 

Now we precede to validation of the correlation relations by comparing our predictions to 
experiments and high fidelity simulations. In 1988 a series of experiments on the natural 
convection in the air (Figure 4-7) has been published in Refs. [14-15], which is considered one 
of the most reliable data on the natural convection to date. 


Validation of correlations. Comparison to experimental data. 


Natural convection of air along a vertical flat plate 



-T. Tsuji and Y. Nagano, “Characten sties of a Turbulent Natural Convection Boundary Layer along a Vertical Flat Plate”, 

Int.J Heat Mass Transfer, 31, pp. 1723-1734 (1988). 

•T. Tsuji and Y. Nagano, “Turbulence Measurements in a Natural Convection Boundary Layer along a Vertical Flat Plate”, 
Int.J Heat Mass Transfer, 31, pp. 2101-2111 (1988). 

Fig.4-7. Turbulent natural convection in the air. Experimental setup. 

A large vertical copper plate was heated to T=60°C, while the surrounding air was held at 
T=15°C. The air was entrained from beneath by the turbulent convection and raised along the 
wall where the hot wire detectors were placed, measuring the temperature and velocity 
distribution. We have simulated the experiment and plotted the results against the experimental 
data (Figure 4-8). 
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3D Simulations of the experiment 



• 3D simulation reproduces well the measured correlations 

• Measurements apparently conflict with Eckert's correlation 


[1] I ‘BUjl and Y. "Otacaetertetics of a Turbulent Natural Convection Boundary Layer a Vertical Flat Hate*, InLJ Heat Mass Tartsfer, 

VdL 31, Ho fi, p*X 1723-1734 (1988). 

[2] I ‘Buji and Y. Nagano> "Turbulence Measurements m a Natural Convection Boundary Layer along a vertical Flat Plate" inti Heat Mass 
Transfer, Vd. 31, No 10, pp. 2101-2111 (1S88). 

[3] E. R. G. Eckert and I W. lacks on, AnalysB of Free-Ccnvectlon Boundary Layer on Flat rtate, National Advisory Committee for Aeronautics 
technical memorandum No. 2207, Washington, September 1950 

Fig.4-8. Turbulent natural convection in the air( High fidelity simulation). 

In addition, comparison to the predictions of Ref. [12] has been performed. The comparison 
is particularly important in view of the fact that correlations proposed in Ref.[12] have been 
widely used in NASA related research for computation of stratification in tanks [18-22]. The 
results (Fig. 4-8) demonstrate a good correspondence of experimental data to our numerical 
simulations, supporting the predicted x' 1/3 scaling of the characteristic velocity, Eq.(29). The 
scaling predictions of Ref. [12] are clearly at odds with these results. 

The next piece of validation comes from simulation of the natural convection of G02 in a 
large cryogenic tank (Figure 4-9). The walls of the tank were kept at T=90K, while the bulk 
temperature of the gas was T=100K. There results of the simulation are plotted on Figure 4-10. 
We see that the predicted x' 1/3 scaling of the characteristic velocity, Eq.(29) is confirmed by the 
results of the simulation. 
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Validation of the correlations. Natural convection of GOx in a tank: 3D 

simulation 
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Fig.4-9. Natural convection of G02 in a cryogenic tank. High fidelity simulation. 


Simulation of natural convection of GOx in a tank. Results 



3D simulations support the predicted scaling of the 
cha racteri stic velocity 

Fig.4-10. Natural convection of G02 in a cryogenic tank. High fidelity simulation 

results. 
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4.2.1. Natural convection in the LH2 tank 


As the further step toward the problem of the temperature stratification in the LH2 tank we 
performed the simulation of the temperature and velocity profiles of the turbulent boundary layer 
flow for the natural convection in LH2 tank (Figure 4-11). The initial (and the bulk) temperature 
of LH2 was taken T„=20.4K. The temperature of the gas-liquid surface (1.4 m from the top of 
the tank) was kept at Ts=23.8K, which corresponds to the partial G02 pressure in the ullage 
PG02 = 2.25atm. The heat flux through the wall was q w =200W/m 2 . The simulation was performed 
for liquid dynamics over 500 sec. At t=50sec after the beginning of the heating, the turbulent 
natural convection flow is expected to have developed. The estimate for the characteristic time t 

of the turbulent flow development is r~ 6 / iu ~ lOsec , where 5~10cm is the boundary layer 

width and u.~1 cm/sec is the scale of turbulent velocity pulsations as measured in the 
simulations. 


Natural convection in LH2 tank. 



boundary layer: apparent width of the turbuler* bounstory layer shrinks 

Fig.4-11. Natural convection in LH2 tank. Velocity and temperature profiles. 

The profile for the temperature in the vicinity of the wall is seen to comply with predictions of 
the theory : T w -T x &0.%K , as calculated using the known heat flux and correlations (A33), the 

width l w &0.5mm as calculated using relation AT <=7^ -7^, Eqs. (31) and (33). Velocity profiles 

corresponding to lines 4 and 5 on Figure 4-11 are seen to collapse on a single graph after the 
best-fit rescaling of velocities by x 1/3 and length by x 8/1 °. The latter scaling deviates a bit from the 
predicted linear dependence of the turbulent outer region width 5, which we attribute the finite 
Reynolds number effect. However the velocity profiles at lines 1,2 and 3 are seen to deviate 
substantially from the prediction of the theory, namely, the apparent turbulent outer region width 
5 is found to saturate at certain value. This deviation could in principle lead to heating the outer 
layer above the temperature predicted by the correlations (29) and therefore should be 
examined carefully. The simulation results presented in Figure 4-12, give the clue to the 
solution. The structure of the turbulent flow in the upper part of the tank is much more complex 
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than in the asymptotic theory of natural convection, which neglects effect of the boundary - the 
gas-liquid interface. Due to this effect, at sufficiently large time t, the descending flow starts to 
develop, which interact with the ascending flow due to the natural convection. We expect that 
this interaction leads to development of the turbulent shear layer which projects deeply in the 
bulk of the liquid far beyond the point, where the velocity of the raising flow vanishes. As a 
consequence, we expect that the actual width of the turbulent flow is much wider than an 
apparent width in Figure 4-1 1 as assessed from the raising flow velocity profile 


Natural convection in LH2 tank. 



Fig.4-12. Natural convection in LH2 tank. Effects of the liquid-gas interface. 

To test our conjecture we numerically computed the distribution of the turbulent kinetic 
energy near the wall for lines 2 and 3 at t=500sec (Figure 4-13). The results of the simulations 
show that in fact the turbulence is spread far beyond the zero-velocity point, leading to an 
efficient heat transport away from the wall. Even if the scaling for this dynamics deviates from 
the prediction of the asymptotic theory of natural convection, the main physical rational which 
lies in the basis of (i) the constant heat transfer coefficient law (33) and (ii) prediction that the 
heating of the liquid is essentially confined to the conductance sublayer of width Iw, is still valid. 
That is, the turbulence efficiently diffuses the absorbed heat in the large volume of liquid, so that 
the main drop of the temperature takes place in the conductance inner region. 

The results of simulation support the foregoing arguments (Figure 4-14): the temperature at 
the wall, taken at t=500sec is found to stay almost constant along the wall and correspond very 
accurately to the predicted value T w —T x »0.SK, the characteristic width of the heated 
sublayer l w =0.5mm corresponds well to the theoretical estimate, too. 
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Natural convection in LH2 tank. 



Fig. 4-1 3. 



"Apparent" turbulent boundary layer width does not correspond to the region of 
turbulent flow, which is much larger due to the interaction with the descenc^ng 
How. 

The turbulent heat transport continues into the region of descending flow 
leading to efficient diffusion of heat in still larger turbulent region. 

Turbulent pulsations vs. average velocity; t=500sec. 
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The turbulent diffusion of heat proves efficient in the Shuttle LH2 tank: the 
temperatureof the wad stays essentially constant as predicted by the theory. 

Ror larger (sis) tank the turbulence is expected to be even lager, leading to 
eve n more efficient dffusion of heat and dec rease of the te mperature in the 
turbulent boundary layer. 


Fig. 4-14. inner sublayer temperature distribution; t=500sec. 


4.3 Temperature stratification in the LH2 tank 


The foregoing results have important implications for the temperature stratification in the 
tank. In fact, the ullage pressure is p u// *33psi=2.25atm during the flight which corresponds to 
the boiling temperature T b0 n a 24K. On the other hand, the predicted near-wall temperature 
T=21.2K corresponds to the pressure of saturated vapor p s , vap =1.1atm. The temperature of the 
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turbulent outer layer is much lower and corresponds essentially to still lower pressure 
Ps.va P =1atm. Therefore, the heating of the liquid due to the natural convection does not create a 
hazard: the Net Positive Suction Pressure condition is satisfied during the flight due to large 
ullage pressure p 0 „ a 33psi=2.25atm. In fact the Net Positive Suction Pressure (NPSP) 
condition for the LH2 tank reads: 

Pa, + Pw - &P/nc ~ Ps, vap (T L (t, x)) > 5 psi (33) 

The hydrostatic pressure p hea d is negligible for LH2. The frictional pressure drop Ap fric is 
determined by drop the thinnest segment of the feedline of the RS-25D engine and is also small 
compared to p u n. We note that since it is a function of the engine only, it is expected to be the 
same for SLS. (Estimation of the frictional pressure drop in tubes is carried out in Appendix.) 

Thus, NPSP condition will hold for the whole LH2 volume of the rescaled rocket perhaps 
only with the exception of very thin layer near the interface of thickness <1.5cm and with 
temperature 23KsT< 24K. This conclusion hold for higher values of the environment heat 
density flow so long as the near wall temperature is (much) smaller than 23. 8K, which 
corresponds to the wall temperature for the critical value of the environmental heat flux 
q .. =950W/m 2 

Te/n’// 


T-stratification in the LH2 tank in the Shuttle 
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Fig.4-15. Temperature stratification in the LH2 tank in the Shuttle; t=500sec. 

The temperature stratification pattern observed in the simulation (Fig. 4-15) corresponds 
well to the theoretical predictions. The warmest layer adjacent to the gas-liquid interface has 
temperature decreasing from T=23.8K to about T=21.2K, as predicted, over the legthscale of 
~1cm, which corresponds to the thermodiffusion length scale l D L=1.5cm calculated for t=500sec. 
The much wider weakly warmed layer at 20.4<T<21 .2 is associated with the contribution of the 
outer turbulent sublayer of the natural convection. Based on the typical velocity in the 
ascending flow and its width (Fig. 4-11) the width of this weakly warmed layer can be estimated 
to be > 8m. To predict the temperature profile in this layer one has to acquire a better 
understanding of the transient characteristics of the turbulent flow pattern of ascending and 
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descending flows near the wall. However, for the scaling predictions for larger tanks it is 
important that the upper limit on the temperature in this layer is expected to stay at or below 
T=21.2K: this is due to even more developed turbulence in larger tanks, leading to stronger 
mixing and diffusion of the absorbed heat in larger bulk of liquid. This latter prediction was 
verified in a simulation of LH2 dynamics in a taller (by the factor 4/3) tank (Figure 4-16), which 
shows that the wall temperature T«21 K is independent of the size of the tank. 


The neaf'wa!! temperature distribution in a large LH2 tank 



The wall temperature is independent of the size of tank as predicted! 


Fig.4-16. Temperature distribution along the wall in a large (40m high) LH2 tank vs. 
the Shuttle LH2 tank (30m high) ; t=50sec, 500sec. 


Our predictions have been verified for the real Shuttle data (Fig. 4-17) for the temperature of the 
incoming fluid in the engine. 


Temperature stratification data in the Shuttle LH2 tank 


LH temperature at 17lnch manifold (data from the flight 
STS 115). Temperature of the incoming fuel at the end of the 
flight rises to T=21.4K. Estimates gives T=2T2forthe heat flux 
q=200W/m A 2 % 
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Fig. 4-1 7 Temperature stratification data vs. predictions for the LH2 tank in the Shuttle. 
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One can see that both prediction of the near wall temperature T=21.2K and the width of the 
weakly warmed liquid layer >8m, correspond well to the data. In fact, the temperature of the 
incoming fluid rises to 21.4K in the end of the flight (the uppermost 1.5cm layer of the liquid at 
T=23.8K is shut-off), and taking the width of the warm layer 10m, we predict the beginning of 
the warm layer inflow at t=400sec. 

In addition, we performed simulation and theoretical estimates of the stratification in Saturn 
LH2 tank, which we compared to the experimental data (Fig. 4-18), as published in [23]. 


Temperature stratification in Saturn LH2 tank 



Fig.4-18 Temperature stratification data vs. predictions for the LH2 tank in Saturn. 

To summarize the main results of our work on the temperature stratification: 

1 . Analysis of the turbulent natural convection correlation relations and the 
condensation-evaporation processes on the gas/liquid interface leads to a fairly accurate 
estimates of the liquid stratification pattern in LH2 tanks in the Space Shuttle and Saturn 
and predictions for larger tanks. 

2. The main risk from the NPSP condition perspective is due to condensation- 
evaporation heating of a thin layer near the gas-liquid interface. 

3. The temperature stratification is essentially invariant under the rescaling of the 
tank. 

4. The interaction of the ascending natural convection flow with the descending 
forced convection flow leads to deviation from the theoretical predictions in the upper 
part of the tank. 

5. However the interaction leads to development of a larger region of turbulence 
where the heat is efficiently diffused. As a consequence the main heating occurs still in 
the conductance (viscous) sublayer where the temperature is independent of the height 
of the tank 

6. The temperature in the turbulent region is of no consequence from the NPSP 
condition standpoint. It is low in the Shuttle and is expected to be still lower in the larger 
tank (SLS) due to stronger turbulence. Numerical simulations of the LH2 dynamics in a 
larger tank support the validity of the theoretical arguments. 
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5. Effect of sloshing to LH2 tank 



Intensive sloshing of cryogenic liquid (LH2) can occur after the start due to vibration of the 
rocket. Two main types of sloshing are shown in Fig. 5-1: (i) the longitudinal sloshing waves 
excited by lateral forces and (ii) the axial sloshing (capillary waves) excited by parametrically 
forced oscillations. 


Fig. 5-1 Two main types of sloshing: (i) Longitudinal sloshing waves excited by lateral 
forces (left); (ii) Axial sloshing (capillary waves) excited parametrically by axial vibrations 
(right) 


5.1. Longitudinal sloshing waves 

To understand main characteristics of the longitudinal sloshing waves (Fig. 5-2) consider a 
rigid tank with ideal inviscid liquid. We suppose that the surface waves are linear and the flow is 
potential 



n ~ 1 mode n- 2rrxxte = 3 mode 

Fig. 5-2 Main modes of longitudinal sloshing 

Let us estimate the frequency of the surface oscillation assuming that the tank has 
cylindrical shape and that the motion is harmonic. In this case the velocity potential ® satisfies 
Eq.(35) 


A<T> = 0 
3 2 <J> c® 

— — + 2 — = 0 at at the free surfacez -hi 2 

dt 2 & oz 


(34) 


! dn = 0 at wall boundaries z--hl 2and r - R,\ 

where h is height of the liquid level and n is a normal to the surface wall, g is a gravity. 
Analytical solution of the eigen-value problem is 

<& mn 0, z, 6) = J x {Kn r / a) cos(m0) (35) 

cosh \XJi!d\ 
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( 36 ) 


which represents eigen-functions and eigen-frequencies 

o£,=^-tanh {XJild) 

a 

where r and 6 are the radial and angular coordinates, a is the tank radius, and X mn 
(m,n=0,1,2..) is a root of the eigen value equation dJ^{X mn r/a)ldr = 0 for r = a, uj is the angular 
frequency . The values of 4 mn for m =1 (cosQ) are A 1n .=^ n . Numerical values of <f„ are: <£, = 1.841, 
12=5.331, = 8.536. . . £, +f -» £, + rr. For m = 1 the potential varies in the angular coordinate as 

cosG, the wave is “up” over half the circumference and “down” over the other half. 

Eigen-frequency of flexural sloshing oscillations is a function of the radius, gravity 
acceleration g, height of the tank h and values A mn . For the first nonsymmetrical mode (m=n=1) 
such frequency is presented as a function of tank radius in Figure 5-3 (left) and the oscillations 
as a function (a,/?) is presented in Fig.5-3(right). It can be seen that for radius R,=4.2m the 
period of the longitudinal antisymmetric sloshing wave is about T osc =3 sec. 

It should be noted that formula (36) doesn't contain density of the liquid, therefore the 
eigen-frequencies will be the same for liquid oxygen tank of the same radius. 


a=r. 



Fig. 5-3. Frequencies of flexural mode oscillations for different tank radios and height 


The longitudinal sloshing is responsible for three main effects: (i) sliding of cold liquid along 
the hot wall; (ii) growth of the area of the liquid surface; and (iii) mixing of the warm surface 
layer with the underlying cold layer. 

(i) At the start of flight the motion of the rocket can induce an longitudinal sloshing wave 
with an amplitude A and period T .The cold liquid in contact with hot tank wall will heat up 

1 wave osc ^ n 


and evaporate. An upper bound on the evaporated mass can be estimated form the following 
condition; the heat accumulated in the wall is spent to heat the liquid and the vapor and to the 
evaporation, Q „ = Q 

1 wall vapor. 



c , (Lau (y))(T al , (y) ~ T s)dy * 


0=interface 


(37) 


Q~, =C,M,XT S -T L )+q L M L +C g M g (T g , vap -T s ),M g =M L . 


35 


Hot wall 


Liquid Wall, T wa „ 



Fig. 5-4. Sloshing brings cold liquid in contact with the hot tank wall 

where we used the fact that the thermodiffusion length in AL 2195 L Q (T osc ) = 7mm is larger than 
the wall width d , ~5 mm. It follows from Eq. (37) that the maximum mass and volume of 
evaporated vapor during t=T osc are 


M. 


evap 


-T s ) 

' <h +C g (T s -T L )+C g M g {T gevap -T s ) 


V = ™r _ «2 y =02m -' 

eavp 9 eavp 

T^tank 


f A ^ 

1 Svavc 

5cm 


for 20K < T evap < 50K. 


(38) 


The pressure jump is calculated using the heat balance equation as follows: 


dE = -pdV + dQ w => ^ M map C p T gas , AV = aAp, 

a = 0.3^ - expansion coeffecient of the LH2 thank 

An — M evap R gas ^ gas 

^ T ullage ,0 + CC rP ullage 

Ap < 0.027 atm for A evap < 1 5 cm, V ullage = 27m 3 , T g - 1 00 K 
A p mm = 0.25 atm for max sloshing amplitude A evap = 1 ,4m 


We see that for the highly improbable hazardous event when the amplitude of longitudinal 
sloshing oscillation at the moment of the engine start (t = - 6sec) equals its maximum value 
A W ave=1.4m the maximum volumes of evaporated vapor and the pressure jump can 
reach V evap = 5.6m 3 and Ap « 0.25 atm for T 0SC =3sec . 

After the start (t > 0) the stationary longitudinal sloshing waves accelerate condensation by 
(ii) increasing the area of the liquid surface (minor effect); (iii) mixing the heated surface liquid 
layer with adjacent cold LH2 with T=T L (major effect). As a consequence, the thickness of the 
heated surface layer decreases but the interface temperature does not change, determined by 
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the partial GOx pressure: T =T s (p GOx ) =23. 8K. The condensation flow rate at the liquid-gas 
interface is limited by the heat balance between the hot gas T g near the surface and the liquid 
surface with T=T : 
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(40) 


The maximum value of the condensed mass M cond induced by the sloshing can be 

estimated assuming complete mixing of the surface layers of cold and warm liquid over each 
period of oscillation. Due to such intensive mixing the width of the heated liquid layer over a 
fractions',’ of the surface area will be very small at end of each oscillation, L l — < >0 (Fig. 5-5). 



Fig. 5-5. Sketch of the temperature distribution before and during the sloshing 


In this case the first term in the right part of Eq. (40) is much larger than the second. 
Therefore, the mass condensed over one period can be estimated as 


M cmd {T 0X )= f S[J cond {t x )dt x *S[{T s -T l ) 
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(41) 


The total condensed mass induced by the longitudinal waves over the time of t is given by 


M:lZ h (t) = M cond (T osc ) 


*S l +S l 

S L 
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4S..+S, j 

v L j 



R. 


(42) 


The total condensation flux is equal to 

max l s,osh — ITtax (0 _ fl nnq kg 

maXJ cond ~ - ^ 

T s (Pghi) = 23.8AT, 24K < T g <50 K near the interface 
It should be compared to the much larger the flux of the injected hot 
GH2: J ghi = = 0.476^ x[Af„„ zz/e /3], N nozde - nozzle nozzle in SLS 


Thus, the longitudinal sloshing waves cannot appreciably change the ullage pressure. 

Effect of sloshing is expected to be even smaller for the SLS since the flow rate of injected 
hot GH2 will increase with larger nozzle number. However, SLS vehicle has higher thrust and 
one can expect stronger vibrations due to adding one or two RS-25 engines. 

Conclusion about the acceleration of condensation by the longitudinal sloshing is confirmed 
by experimental data [23] presented in Fig. 5-6. The pressure in the ullage space increases 
from 15psi to p o =35 psi due to hydrogen pressurization with temperature 40K during 15sec. 
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After 50 sec lateral vibrations are observed with frequencies f=0.95Hz (non-resonance 
influence) and 0.74Hz (resonant vibrations stimulate intensive sloshing). The ullage pressure 
slowly decreases in the case of the vibrations with f=0.95Hz and sharply drops in the case of 
the resonant vibrations with f=0.74Hz. 
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Fig. 5-6. Sketch of the tank and the internal instrumentation (a), pressure response for 
the sloshing test (b), wall temperature (c) and ullage space temperature (d). 


The condensation flow density rate can be estimated from the equations 


p{t)V = p 0 V - j cond S ■ t, p Q V = p 0 VR H 2 T 0 , 
py — ( Po V ~ Jcond^ ' f| ) Rtf 2 ^1 = 4 1 24 ^r), 
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J cond 


T x 

Po — - 

\ y o 


Pi 


R H2 T x S-t { 


It follows from Eq. (44) 

j cond, \ ~ 4 5 • * for resonance (f=a74Hz ) 

for p 0 =35psi. p,=20psi, At=20sec,V ullage =0.35V lank =0.63m \ S=1.6m 2 , 


J cond 7 —3-3' i ^2g gc f° r non-resonance (f=095Hz) 

We use the data of Eq. (45) for estimation of the condensation rate for LH2 tank 


(44) 


(45) 


J cond ~ JcondS ' V -* f* ~0A6 for resonance influence (f, =0.74Hz) 

Jcond ~ JcondS ' facfi)' ~ 0014 A for non-resonance influence (f 2 =0.95Hz) 

The foregoing estimation = 0.025|§for the condensation rate for intensive sloshing 
given by Eq. (42) agrees with the data of Eq. (46). 
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5.2. 


Axially excited sloshing: generation of droplets 


The thrust oscillations initiate axial excitation and lead to Faraday instability of the liquid 
surface. Intensive oscillations of these modes can result in the liquid level discontinuity and 
generation of the liquid droplets under the liquid-gas interface. This process will change heat 
balance on the interface. Evaporation of the cold liquid droplets into the hot ullage gas can 
increase the ullage pressure during the beginning stage of the flight. At the same time enlarging 
the cold surface can accelerate condensation. Since the rate of evaporation of liquid hydrogen 
depends on partial pressure of gaseous H2, in the case where He is used as the pressurization 
gas we expect enhanced evaporation and consecutive increase in a ullage gas pressure. For 
gaseous Hydrogen pressurization these droplets will enhance the condensation rate. 

Theoretical treatment of the droplet ejection is given in many articles. Linear theory is 
described by dispersion relation obtained from equation (34) with the same boundary 
conditions, where for small wavelength capillary waves we should take into account surface 
tension a : 


*> 2 = 


^ , a 
gk+ — k' 

P J 


ia.nh( kh) 


(47) 


where k=2rr/A, A as the wavelength, is the wave number and p is the liquid density. At small 
k«1 this equation is identical to equation (36): 

gA " m -1anh(A nm h/a) (48) 


<y;„, =- 


a 


The nonlinear dynamics of the droplet formation and ejection is extremely complicated and 
strongly depends on the rocket vibrations spectrum. Here we carry out a simplified analysis of 
this dynamics. The typical diameter of ejected droplets is about [24,25] 

djrop - K / P L ) m co~ v \ (0 = Inf (49) 

where p L =70kg/m 3 and o L =0.0017N/m are liquid density and surface tension, co is vibration 
frequency (a=2nf vibr ). Characteristic values of d drop are shown in Fig.5-7 (left). 



Fig. 5-7. The typical diameter of ejected droplets as a function of the thrust 
oscillation frequency (left) ; generation of droplets under action of supercritical vibration 
acceleration £ d >0 [26]. 
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Each projectile area of diameter d proj ~5d drop ejects the droplet with a rate given by [24,25] 

dn/dts i ) = 0. 04f 7*0 [s' l ],Sj = (a viir -a d )/a d , 

a d = 0.26(<r / p) m o m 

Here a, co are vibration acceleration and frequency, a d is the critical driving acceleration for 
the onset of droplet ejection. The droplet ejection occurs when the frequency and amplitude of vibration A v ,br 
satisfy the condition: 

> 4 / or A, hr >4 = °- lmm(\00Hz / / V/A ,) 2 ' 3 , 

A c = 0.25 mm for f = 25 Hz, A c = 65pm for f = 266 Hz 

Sufficiently strong vibrations can occur only after the start. At this point the thermo-diffusion 
length in the liquid has increased to the value of L q (100sec)=7mm during the prepress (t * 

lOOsec). This value is much greater than the typical drop diameters: 
L d — Imn » d ii71p forf >50 Hz. It follows that the inner temperature of generated droplets 

T d op is close to the surface temperature of the liquid T s (p GH2 ). Therefore, the liquid mass 

evaporated due to the vibration is very small. Indeed, the mass evaporated from the surface of 
single droplet s drop during the droplet life-time t drffe can be estimated as 


, (t fW''*'* • f/W/ < ^ c»nd( l presf S dm r _ ffrop 

evap\ l drlife ' ^ drop I Jevap' 1 * 1 1 drlife ’ * drlife \{ 9 

o L V g 
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(52) 


where h drop is the maximum height an ejected droplet rises above the surface (see Fig. 6-7). 
The total evaporated mass is given by 

dNT 1 

l evaop V drlife J 


'±2E- t = m cond (t d ) -— p - t. 


evap 


dt 


pro i 


dt ■*" dt 


(53) 


Here dNf° p / dt is total number of droplets ejected over unit time from the liquid surface S L . 
Taking into account Eq. (52) and condition d proJ ~5d dmB ,s drilD / s Dmi - 4/25 we conclude from 


A drop 9 ° drop ' ° proj 


Eq.(53) that the total evaporation mass is a negligible quantity 


^ evap ~ J evap V press drlife ' 
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L drlife (54) 

J evap * M evap ' ' = 0. 1 6 J evap (T prepre J < 0.8g / SCC 

In the estimation we used the results of calculation presented in Sec. 2 (see Fig. 2-4) 
J (T )<5e/sec . 

An appreciable effect can occur only as a result of combined action of intensive longitudinal 
and axial sloshing vibrations (Fig. 5-5, left)). The maximum value of the condensed mass 
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induced by sloshing waves can be estimated assuming intensive mixing of surface layers of 
cold and warm liquid during each period. Due to such mixing the thickness of some area of the 

heated surface layer S', will very small at end of each oscillation, l_ L — >0 (Fig. 5-5). Droplets 
generated by the axial sloshing will have T near T l , except for a very thin surface layer l_ L — >0 
near the surface having T=T s (Pgh 2 )- In this case the maximum value of the condensed mass 


m cond on the surface s drop °f a single droplet during the droplet life t drffe can be estimated as 
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(55) 


and the total condensation flux from all the ejected droplets is 
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Using S', ~S, /3 and t Mife * f2h lmp 1 g «0.02sec we obtain from Eq.(56) that the total mass 
flow rate of condensed vapor on all ejected droplets is about J'f° p < 0.09 kg t sec . We notice that 

the mass flow rate of injected hot GH2 is much more: J GH2 = 0.35^- = 0.476^ *[N noz2le / 3] . 

As a result of the condensation the ejected droplets will be heated and falling back into the 
liquid will increase the temperature of the liquid surface layer. The width of the heated surface 
liquid layer can be found from the heat balance equation 

•C [V Cw: fc - r, )) ' = C,A (: T,-T l )S l -I**, 

,t*u<d (j) - . j 1 - ^ i - 

1 heated VJ ^ ^ j\ l l drhfe ) fL d,-»p? L Jrhfe ^ ^ > 

^ heated flight ) — 5L D (t flight) = 2 .5 cm for fc, Vt?i < 1 / t ithk , t Jrhfe < 0.02 sec 


(57) 


Thus, combined action of longitudinal and axial sloshing vibrations can lead to formation of 
relatively thick hazard surface layer of heated liquid with temperature T a 23.8K, thickness 

L~7.5cm, volume , = S i:mk /!;f irJ / 2 = 2m' and mass = V^p, = 140% . Total mass of 

LH2 in the tank is 105,000kg 

SLS vehicle has higher thrust and expected vibrations due to adding one or two RS-25 
engines. At the same time the mass flow rate of injected hot GH2 increases with the engine 
number. All the effects under consideration are taking place in a narrow layer adjacent to the 
interface. Therefore the influence of these effects on both the mass and the heat balance in the 
tank will not be critical. However, the effect of heating of relatively thin layer to T«23.8K can 
arise in LH2 SLS tank because possible intensive longitudinal and axial sloshing vibrations . 

6 Conclusion 


We performed theoretical investigation of main physics processes in LH2 tank during pre- 
starting prepress and rocket flight that can lead to possible hazards in the case of rescaling of 
main parameters of Shuttle with regard to proposed LH2 tank designs for SLS with 5 engines 
(the situation with 4 engines is less critical). Our investigations based both on theoretical 
analysis and on MNM and 3D simulations allow us to make the following conclusions 
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1 . Hazard analysis for LH2 ET shows that the pressure control bandwidth remains almost 
the same for the SLS. 

2. The MNM can reproduce accurately the pressure during pre-press and blow down. 

3. For the GHe and GH2 flow rates and the initial ullage volume rescaled by 5/3 for the 
SLS LH2 ET, no new hazards are found for the tank performance; the temperature stratification 
is more pronounced in the rescaled SLS tank. 

4. The bubble formation at the tank wall does not scale with the tank length but it's effect is 
still relatively small and will not considerably affect the performance. 

5. Effect of condensation and evaporation are also relatively small. 

6. The temperature stratification of the liquid reveals a structure consisting of three layers: 

(i) The upper hearted layer of thickness less than 0.6in=1.5cm at the end of the flight and 
with temperature 21.2K<T< 23. 8K; 

(ii) The intermediate warm layer of width <8in=20cm at the end of the flight and with 
temperature 20.6K<T < 21.2K; 

(iii) The lower warm layer with 20.4K<T<21.6K of width ~10m at the end of the flight. 

In (i), both the temperature and the width of this layer are invariant under the 5/3 rescaling of 
the tank length. In (ii-iii), the sub-layer is formed due to the rising warm liquid flow with T=21.2K 
forming in a thin boundary layer near the wall and accumulating near the interface. These sub- 
layers will increase by a factor of 5/3 under the scaling but its temperature will be T< 21. 2K. 

7. The NPSP condition holds for the rescaled rocket with the exception of the upper 
hearted layer of thickness less than 0.6in=1.5cm . Maximum mass and energy of such 

overheated liquid with T=23.8 is less than 30kg and 2.5x10 J. Total LH2 mass and energy in the 

8 

filled tank are 105270kg and 3.4x10 J. Therefore, it does not play the role in the ullage 
temperature stratification and pressure dynamics. 

8. Intensive longitudinal sloshing waves induced by the laterally excited rocket vibration 
change weakly the temperature liquid and ullage stratification and the ullage pressure during the 
flight. Maximum rate of condensation flow arising due to the strong mixing of the warm surface 

liquid layer with the cold deep layers is limited by value of J\ = 0.023^- . The mass flow rate 

of injected hot gas (GH2 or GOx) is much more JcaXO = 0A76-& x(n m2de /3) , n nozz i e is nozzle 
number in the SLS. 

9. The vertical vibrations with the amplitude A vibr and the frequency f vibr such as 
Aib- > A -0.\mm(\00Hz / f vlbr f n induce generation of the warm droplets under the liquid-gas 

interface. Effect of axially excited vibration sloshing on LOx tank is extremely small as the 
temperature of generated warm droplets is very close to the interface temperature and 
condensation is depressed. 

10. Combined action of intensive longitudinal and axial sloshing vibrations can induce 
relatively great condensation and formation of relatively thick hazard surface layer of heated 
liquid with temperature T=T s (p G ox)- The rate of the condensation can 

reach J'^ p = 0.09kg /sec and thickness, volume, and mass of surface heated layer can reach 

L L =7.5cm, V=4m 3 , and 280kg. Total LH2 mass and energy in the filled tank are 105270kg. 

11. The longitudinal and axial sloshing waves cannot appreciably change the ullage 
pressure. Effect of sloshing is expected to be even smaller for the SLS since the flow rate of 
injected hot GH2 will increase with larger nozzle number. However, SLS vehicle has higher 
thrust and one can expect stronger vibrations due to adding one or two RS-25 engines. 
Combined action of intensive longitudinal and axial sloshing vibrations can lead to formation of 
relatively thick hazard surface layer of heated liquid with temperature T=23.8K, thickness 
L=7.5cm and volume (mass) V=4m 3 .(280kg). 

12. SLS vehicle has higher thrust and expected vibrations due to adding one or two RS-25 
engines. While the mass flow rate of injected hot gases increases with the engine number. All 
the effects under consideration are taking place in a narrow layer adjacent to the interface and 
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the influence of these effects on both the mass and the heat balance in the tank will not be 
critical. However, the effect of heating of relatively thin layer to T=23.8K can arise in LH2 SLS 
tank because of possible intensive longitudinal and axial sloshing vibrations. 


7 Appendix 


The pressure drop in a tube with turbulent flow due to the wall friction can be estimated as 
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(58) 


Re > 4 x 1 0 7 for J L > 200kg / sec and r luhe > 0.2 m (turbulent condition) 

Requirement for the feedline follows form Eq.(58) 

p fnc < 0. 1 atm for J=200kg/sec (3engines) at tube radius r lube >7.5cm(3in) 
or for J=333kg/sec (5engines) at tube radius r tube >9cm(3.5in) 


(59) 


Thus, in order to keep the same pressure drop in the feedlines in scalable SLS LH2 tank it is 
necessary to increase the tube radius on 20% or 10% for nozzle number N =5(4): 
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= 1 .2(1 . 1) for n nmzle = 5(4) and U = l luhefl . 
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